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(Preliminary  draft) 


W.  Kohn,  J.  H.  Jensen,  and  P.  Chang 
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ABSTRACT:  We  consider  an  atom  scattered  inelastically  from  a  surface  and 
derive  a  formula  for  the  exact  semiclassical  limit  of  the  expectation  value  of  an 
arbitrary,  smooth  function  of  the  scattered  atom’s  final  momentum.  The  formula 
expresses  this  semiclassical  limit  in  terms  of  equilibrium  correlations  for  the  atoms 
which  compose  the  surface  and  quantities  which  can  be  calculated  from  classical 
scattering  trajectories.  We  give  numerical  results  for  a  simple  one-dimensional  ex¬ 
ample  and  compare  these  results  with  those  given  by  the  trajectory  approximation. 


I.  Introduction 


Atom-surface  scattering  can  be  used  to  obtain  information  about  the  structure  and 
dynamics  of  a  surface  and  about  the  atom-surface  interaction  potential.1  When  the  initial 
energy  of  an  incident  atom  (hereafter  called  the  “incidon”)  is  large  and  its  de  Broglie 
wavelength  is  small,  the  semiclassical  limit  provides  an  accurate  approximation.  In  this 
paper,  we  derive  an  expression  for  the  exact  semiclassical  limit  of  expectation  values  of 
functions  of  the  incidon ’s  final  momentum  in  terms  of  classical  trajectories  and  equilibrium 
correlations  for  the  atoms  composing  the  surface.  Since  there  are  well-established  methods 
for  calculating  classical  trajectories  for  realistic  surfaces,2  the  formula  provides  a  practical 
method  to  obtain  the  semiclassical  approximation  for  quantities  observed  in  scattering 
experiments. 

We  illustrate  the  theory  with  a  one-dimensional  model  consisting  of  an  incidon  scat¬ 
tering  off  a  harmonic  oscillator,  giving  numerical  results  for  the  semiclassical  momentum 
shift  and  momentum  uncertainty.  We  compare  these  results  to  those  produced  by  the 
trajectory  approximation  which  is  frequently  used  in  the  study  of  inelastic  atom-surface 
scattering.3-9  We  find  that  the  trajectory  approximation  does  not,  in  general,  give  either 
the  correct  semiclassical  shift  or  uncertainty.  However,  if  the  incidon-oscillator  interaction 
is  adiabatic  or  if  the  maximum  displacement  of  the  harmonic  oscillator  from  its  equilibrium 
position  is  much  smaller  than  the  range  of  the  interaction  potential,  then  the  trajectory 
approximation  accurately  predicts  the  momentum  uncertainty. 

II.  Definition  of  Semiclassical  Limit 

Many  approaches  to  surface  scattering  have  gone  under  the  name  “semiclassical.”  Of¬ 
ten  the  term  means  that  the  approach  employs  both  classical  and  quantum  mechanical 


»*»  IViViS  ».*  ».t  , 


methods.  Alternatively,  semiclassical  can  imply  an  expansion  in  powers  of  h.  It  is  in  this 


latter,  more  precise  sense  that  we  use  the  word. 


We  consider  a  d-dimensional  quantum  system  consisting  of  an  incidon  and  a  surface 


substrate  of  harmonic  oscillators.  The  Hamiltonian  is 


H  =  +  ^(x»  6)  +  9  £(*?  +  w? £?)> 


where  p  is  the  incidon ’s  momentum,  x  is  the  incidon ’s  position,  m  is  the  incidon ’s  mass,  and 


V  (x,  £,)  is  the  atom-surface  interaction  potential  which  is  taken  to  be  a  smooth  function 
of  x  and  o{  and  £,•  are  coordinates  for  the  normal  modes  of  the  surface  atoms,  obeying 
the  commutation  relation  aj\  =  ihSij ,  and  w,-  is  a  normal  mode  frequency.  We  assume 
that  V (x,  £,)  — ►  +oo  as  x  ■  n  — ►  — oo,  and  V(x,  &)  — ►  0  as  x  •  n  — »  +oo,  where  n  is  a  unit 


vector  perpendicular  to  the  surface. 


At  an  initial  time,  tQ ,  the  surface  is  in  thermal  equilibrium  at  a  temperature,  T  —  hr, 


and  the  incidon  is  described  by  a  wave-packet  of  the  form 


V>(x)  =  (271A7)  *  exp 


i _ [  (x  ~  x0)2  po  ■  x 


where  Xq  is  the  incidon’s  initial  mean  position  and  po  is  its  initial  mean  momentum, 


satisfying  po  •  n  <  0. 


At  a  later  time,  f,  the  expectation  value  of  a  smooth  function,  Q( p),  is  given  by 


<  Q(P)  >  (0  =  Q{t  ~  to,XQ,Po,~1,T,h). 


IS 


— — 1 

Since  the  impact  parameter  for  the  collision  is  generally  unknown,  we  average  <  Q(p)  >  if  ’v, 

H  8 

over  the  components  of  Xq  which  lie  parallel  to  the  plane  of  the  surface.  Indicating  this 


average  with  a  subscript  a,  we  have 


<  Q( p)  >a  (*)  =  Qa(t  -to,X0  •  n,p0,7,  r,  h). 


(4) 


W 


The  final  classical  value  for  Q  is  defined 

Qfd Po.O  = 


lim  (  lim  (  lim  (  lim iQa(t  ~  *o»xO  ‘  n.Po,'/,  r,  ft) ) )  )> 

>00  \x0-n— ‘oo '(i— tn1)— >oo 'A— '0 


(5) 


7—oo  \x0-n—  oc^(t-to) 
while  the  semiclassical  shift  for  Q  is 
AQ(po,r,fc)  = 


ft  lim  f  lim  f  lim  (  lim  — Qa(t  -  to»x0  •  n,p0,7,r,ft) ) ) ) . 

7— oo\x0-n— oo \t— to)— *00  'fi— 0  on  11 J 

Finally,  we  define  the  semiclassical  limit  or  approximation  for  Q  as  Qsc  =  Q fc  +  A Q. 


(6) 


If  the  incidon’s  initial  energy  is  large  compared  to  the  energy  of  a  typical  phonon  and 
the  incidon’s  de  Broglie  wavelength  is  small  compared  to  the  distance  scale  over  which  the 
interaction  potential  varies  substantially,  then  we  expect  QSc  to  be  a  good  approximation 
for  the  final  expectation  value  of  Q.  Since  we  have  taken  the  surface  temperature  to  be  of 
order  ft,  the  initial  energy  must  also  be  large  compared  to  k^T  with  kg  being  Boltzman’s 
constant. 


III.  Derivation  of  Formula  for  Semiclassical  Limit 


To  obtain  an  expression  for  the  semiclassical  limit,  we  use  the  Wigner  distribution 
function10  defined  by 

ddaDa; 


/(*,p 


exp 


~l(P'a  +  SaiQ!*) 


N 

c 

i=l 


(7) 


(27rft)(d+iV) 

/  1  1  A  1  A  1  \ 

x  p  -t-  -a,x  —  -a,  &  +  -a;,  £,■  —  -ai,  t  J  , 

where  Dai  =  and  p(x,x',  &,  £'-,i)  =<  x,  ^^(Olx',  >  is  the  density  matrix. 

The  expectation  value  of  <?(p)  can  be  written 


<  Q{ P)  >  (0  =  J  ddxddpD(iD<7iQ{p)f{x,p,(i,ai,t). 


(8) 


The  equation  of  motion  for  /  is 

/  d  pi  d  d  o  *  d  \  .  .  . 

(*+m 5* ' + ' «  Wi  - “Ui  ad  /(x' P' {il  "•  0  = 

-Y^j(x,zi,ii,\i)^ Tl[jdpU))  [n(2 da(k))  j/(x’p,ei',a,,)’ 

where  the  sum  is  over  all  positive,  odd  L,  with  L  =  J2j=l  lj  +  Hk= i  ^ k  where  lj  and 
A h  are  nonnegative  integers;  repeated  indices  are  summed  unless  they  are  contained  in 
parentheses,  and  the  coefficients,  J,  are  given  by 


The  initial  condition  is 


/(x,p,&,<Ti,t0)  = 

fsur(Zi,0i) 

(„f,\d  eXp 


(x  -x0)2  2^(p  -  po)5 


where 


/,Ur(f, exp[-tanh(^)  (tt* + Sr)]-  (12) 


The  classical  trajectory  is  defined  by  the  equations 
pc/(i)  =mxcl(t) 

•  cl(f\  _  _  dV  (*>  &)  I 

P  '  x=xe,(0;f«=f*'(0 

,  .  ,  (13) 

of(t)  =«'«) 

i cl  It) _ 2  cel  ft)  _  &)  i 

°i  (*)  I x=xe'(0;£, -#'(<)’ 

and  the  initial  conditions  xc,(fo)  =  xo,  pc,(*o)  =  P0>  £,-*(*o)  =  °>  and  We 

assume  the  incidon  does  not  stick  to  the  surface  so  that  xcl(t)  •  n  — ►  oo  as  t  — *■  oo. 

Defining  a  distribution  function,  g,  so  that 

ff(r,s,ei,  ?,-,«)  =  (ft)<f+iV/(x,p,&,<7j,t),  (14) 
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i’,: 
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li 
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s» 
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ffi 
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O' 
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where 


equation  (9)  may  be  rewritten 
d  S{  d  d  o  d 


r=^=[x-xc/(t)l 
s=-L(p-p''(1)i 
c,  =^!fi  -  «?(«)! 


f  d  S{  d  6  2  d  1  dV ,  d  1  dV ,  d  \  ,  ^ 

{di+mfrri+qidei  ~ Ui€idsi  +  y/KdxMdsi  +  y/kd^cld(Xi)9{T,8,€i,<;i,  )  ~ 


(15) 


(16) 


where  the  symbol  “|c/”  indirates  that  a  quantity  is  evaluated  for  the  classical  trajectory. 
The  initial  condition  for  g  is 

g(r,s,Ci,Ci,t0)  = 

_*.»].  (17> 

where 

=  n  ;  tMlh(^)  «!>[-  tanh(^)  ("Ww  +  ^)]  ■ -  (18> 

and  the  coefficients,  J,  can  be  written 

r  d  «  .  i.  -i  r  N 


Jfe,  ^i)  ~ 

OO  r  d 


5  \5(n) 


(19) 


V’tx.fiJIrf, 


where  A  =  bm  +  5Z^=i  0n,  with  bm  and  9n  being  nonnegative  integers. 

From  equations  (16)-(19),  one  can  show  that  g  has  an  expansion 

OO 

S(r,S,6i,?j,0  =  ^(ft)2ffn(r,S,£j,Cj,0, 


(20) 


where  gn  is  independent  of  h  and  has  the  properties 


J  ddrddsDtiD<;igQ{r,  s>  et>  ?i.  0  =  1 

J  ddrddsDf.i D<:ign (r,  s, <?i,  *)  =  0,  n  /  0  (21) 

5„(-r,-s, *)  =  (-l)nff„(r,s,e,-,  ?,-,*)• 

The  equations  of  motion  for  the  gn  are  found  by  substituting  the  expansion  for  g  from 
equation  (20)  into  equation  (16)  and  collecting  terms  of  the  same  order  in  h.  Using 
equations  (8),  (14),  (15),  (20),  and  (21),  the  ft-expansion  for  <  Q[ p)  >  can  be  shown  to 
be 


Q{ p)  >  («)  =  ^(p)lc<  +  ^[-^7,le/  J  ddrdd3DeiDgigi{rt8,eit(:i,t)3j 


+ 


\ddjrdpk  'c/  /  d<lrddsDeiD<:i9 +  ®(ft2). 


(22) 


Now  consider  a  function,  /,  which  satisfies 


ffL+PlA 

\dt  mdxi  ld£i 
dV(x,Zi)  d 
dxi  dpi 


)/(x,p,ei,(Ti,t)  = 

dv(x,  ii)  d  i?,_.  _  ,  _ 

,  j  /(^iPi  £»»0»*O» 


+ 


d£i  da{ 


(23) 


with  the  initial  condition  /(x,  p,  <r,-,  to)  =  /(x,p,  £,-,<7,-,to)>  Equation  (23)  is  the  clas¬ 
sical  transport  equation  for  a  phase  space  probaoility  distribution  corresponding  to  the 
Hamiltonian  H  and  is  identical  to  equation  (9)  except  that  the  terms  with  L  ^  1  do  not 
appear. 


The  function  g,  defined 


sKr,s,€j,<Ti,t)  =  (ft)d+^/(x,p,&,Oj,t), 


(24) 


satisfies  an  equation  identical  to  equation  (16)  but  without  the  I  ^  1  terms  and  has  an 
ft-expansion 


CO 


g(r,s,  =  ]>^(ft)  2  £n(r,s,  £;,?,•,*). 

n=0 


(25) 


g  also  has  the  properties  given  in  equation  (21). 


The  equations  of  motion  for  go  and  go  are  identical,  and  hence  g q  =  g q  for  all  times. 
Although  the  equations  for  g\  and  31  are  different,  it  can  be  shown  (see  Appendix  A)  that 

J  d?rddsDciD<;igi(r,9,ei,<;i,t)sj  =  J  ddrddsDeiD<;igi[r,8,ei,<;i,t)sj.  (26) 

From  equation  (22),  we  then  have 

<  <3(p)  >  (0  =Q(p)|c,  +  ^*|gl|ci  j  ddrddsDeiD<:igi{Tta,ei,<;i,t)sj 

+  J  dV’<isI>elI><:i3o(r,s,ei,?l-,Oaj5jk]  +  0[h?)  (27) 

=  J  ddxddpDZiDoiQ(v)f(x,-P,Zi,oi,t)  +  0(h2). 

Therefore,  the  classical  nhase  space  distribution,  /,  gives  the  correct  quantum  expectation 
value  for  Q  up  to  order  h. 

We  define  a  set  of  classical  trajectories  by 

p<r(t,x,,p/,  £j,<rj-)  =mxlr(t,x,1  p7,  £j-,aj) 


^r(f,x',p',^,<r-)  =  - 


ay(x,&), 


o|r(t,X,,p',^,CT-)  =€ir(t,x',p',^,0-) 

(7jr(t,X,,p,,eJ,<T|)  =-W^)^)(t,x',p,,e>')  - 

where  the  derivatives  of  V^(x,  &)  are  evaluated  at 

x  =x*r(f,x/,p/,  ^|,o() 


dV(x ,  e,) 


and  with  the  initial  conditions 


6  =£-r(f,x',p',£,<r,'). 


X<r(t0,X,,p',e-,<T')  =X' 
p‘r(to,x\p\&^)  =p' 
£!r(*o>x'>p'>^X)  =€i 

<7ir(t0»x'.p'.^-.^)  =°'i ■ 


w' 
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S!>; 

Si! 
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The  solution  to  equation  (23)  can  then  be  written 

/(x,p,ei,<Ti,i)  =  f  ddxUdp'Dt'Mnxl,v,,M,t0) 

x  <5d[x  -  xtr(t,x',p',  ~  P<r(f>x/>P,»  (31) 

x  fl  {6^  -  ^’x'’p 

i=l  1 

which  implies 

<  Q(p)  >  (0  =  I  ddx'ddp' D^Da[Q{ptr f[x\ t0) +  0{h2).  (32) 

Expanding  about  the  initial  values  for  the  classical  trajectory  of  equation  (13)  yields 

~clt 


<  Q{ p)  >  (0  =  O[pc/(01  +  \~^'[cl8‘lx'i  +  2  dp f  'dfi2pi 


+  3@'^  +  l@^  +  01^ 


(33) 


where 


lx'=Xo 


a2<?  a2ofptr(t,x'>P0, o,o)) 

axf,c'“  a*? 

82Q,  82Q[pfr(t,Xo,p;, 0,0)1  ■ 

S^ld  - - Stf lp  =p» 

d2Q.  a2Q[pir(<,x0,P0.^.°)l| 

5{Slci  =  l5:=0 

d2Q,  a2Qlp<r(t,Xo,Po,0,a()]( 

l*=0. 


(34) 


and 


(35) 


62x\  =  [  d?x'd*p'D?iDo,if{x',p\&(7i,to)[xi  -  (x0 ) j]2  =  hi 

S2p[  =  J  dJz'dJp'D('iDo'J(x',p',('i,<r'i,to)lp'i  ~  (Po)if  = 

S2('i  =  f  =  2^coth2%r 

<Vi  =j  d^^p'D^DcUV-P'’  «1X. ‘0)«?  =  ^rco*h^7- 

In  order  to  obtain  the  semiclassical  shift,  AQ,  we  must  average  equation  (33)  over  the 
impact  parameter,  differentiate  with  respect  to  S,  and  take  the  limits  indicated  by  equation 


1 


a 


K  jj 

yj 
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(6).  Doing  this  causes  the  first  three  terms  on  the  right  side  of  equation  (33)  to  vanish,  as 
well  as  those  of  order  A2  or  higher.  The  first  term  vanishes  since  it  is  independent  of  h ,  and 
the  third  term  vanishes  since  8^p\  goes  to  zero  as  7  tends  to  infinity.  The  components  of  the 
second  term  perpendicular  to  n  vanish  upon  averaging  over  the  impact  parameter,  while 
the  component  parallel  to  n  vanishes  since  the  final  momentum  of  the  incidon  becomes 
independent  of  Xq  •  n  as  (t  —  t^)  —►00  and  Xq  •  n  — *■  00.  Thus,  the  semiclassical  shift  is 


given  by 


where 


ld2<?a,  c2ci,\^9±\  S2  J 

A<?  ~  2  d£'?  \f°6  ^  +  2  do'-2  fc6  *’ 


vv  ^ 

%|/c  =  lim  (  am  (!^u)  ) 
xo-n— 00  00 '3^2  '  a' 

d2Qa  _  (  (S^Q_  \  \  (37) 

with  the  subscript  a  indicating  the  average  over  the  impact  parameter  and  the  subscript 
fc  indicating  a  quantity  evaluated  for  the  final,  classical  scattered  state  after  the  limits 
xq  •  n  — ►  00  and  (t  —  (q)  — *■  00  have  been  taken. 

Equation  (36)  may  be  rewritten  in  terms  of  coordinates.,  y,-,  which  represent  the  dis¬ 
placements  of  the  surface  atoms  from  their  equilibrium  positions,  as 

AQ  =  <  ViVi  >0  +i4%'/c  <  Ms  >0’  (38) 

where  the  correlations  for  the  surface  atom  coordinates  are  evaluated  in  thermal  equilib¬ 
rium.  Equation  (38)  has  a  simple  physical  interpretation.  Consider  the  classical  scattering 
of  an  incidon  off  a  surface  where  the  surface  atom  coordinates  initially  obey  a  Gaussian 
probability  distribution  with  the  property  that  it  implies  atom-atom  correlations  identical 
to  the  equilibrium  quantum  mechanical  atom-atom  correlations.  Then  the  average  value 
of  Q(p)  for  the  final  scattered  state  will  be  equal  to  the  expectation  value  of  Q(p)  for  the 
corresponding  quantum  problem  to  order  ft.  Approximations  using  classical  trajectories 
with  quantum  initial  conditions  have  been  applied  to  molecular  scattering11  and  to  the 


calculation  of  relaxation  rates  for  particles  interacting  with  heat  baths.12,13  Equation  (38) 
justifies  this  approach  for  atom-surface  scattering  in  the  semiclassical  limit. 


We  finally  note  that  if  2k qT  hu>{  for  all  w,-,  then  AQ  a  T  for  any  Q. 


IV.  One-dimensional  Example  with 


Comparison  to  Trajectory  Approximation 


As  an  application  of  the  results  of  the  previous  section,  we  consider  a  one-dimensional 
example  where  the  surface  consists  of  a  single  harmonic  oscillator.  The  Hamiltonian  is 

=  ^  +  ^  +  (39> 
where  q  =  \[Mo  and  y  =  £/VM  with  a  and  £  being  the  oscillator’s  normal  coordinates, 
and  M  is  the  oscillator’s  mass.  Equation  (38)  gives 

1  d2p ,  0  1  d2p ,  .o 

A?  =  2^/2  l/c  <  y  >0  I fc  <  V  >0>  (4°) 


2  3y'2 


6'2P  =  <  y 2  >0  +(^|/c)2  <  ^  >0’  (41) 

where  S2p  =<  (p  —  pcl)2  > f  +0(ft2)  is  the  momentum  uncertainty  of  the  final  incidon 


state.  Equation  (38)  also  implies 


.  ~  dQ(p)  |  ,2 

AQ  ~^r'/cAp  p 


^  /.)  W  («) 

where  <52Q  =<  [Q(p)  —  Qc/]2  >/  +0(ft2).  Equations  (42)  and  (43)  imply  that  the  semiclas¬ 
sical  shift  and  uncertainty  for  any  quantity  may  be  found  from  the  shift  and  uncertainty 
for  p  and  the  final  value  of  pcl. 


The  trajectory  approximation  has  been  developed  to  calculate  the  final  energy  distri¬ 
bution  of  the  incidon  (and  for  a  flat  surface  the  momentum  distribution  as  well).  Although 
predictions  based  on  the  trajectory  approximation  have  been  compared  to  experiment  with 
some  apparent  success,^-®  a  compelling  theoretical  justification  for  it  has  not  been  given. 
To  test  its  validity,  we  compare  the  semiclassical  momentum  shift  and  uncertainty  pre¬ 
dicted  by  the  trajectory  approximation  to  the  exact  results  calculated  from  equations  (40) 
and  (41). 

The  essence  of  the  trajectory  approximation  is  that  the  incidon  is  assumed  to  follow  the 
classical  trajectory,  while  the  surface  atoms  are  treated  quantum  mechanically  with  a  time- 
dependent  interaction  potential,  V(xc/(t),  £,].  The  final  energy  distribution,  PS{E ),  for  the 
surface  is  then  calculated,  and  the  final  energy  distribution  for  the  incidon  is  inferred  from 


the  relationship 


P,nd(E)  =  J  dE'P',{E')P,(E0  +  E’~  E ), 


which  follows  from  energy  conservation,  where  P'S{E')  is  the  initial  energy  distribution  for 
the  surface  and  Eq  is  the  initial  incidon  energy. 

Using  arguments  similar  to  those  of  section  III,  we  can  obtain  expressions  for  the  semi¬ 
classical  limit  of  the  final  energy  shift  and  energy  width  in  the  trajectory  approximation. 
For  simplicity,  we  confine  ourselves  to  the  particular  system  described  by  We  define 


a  set  of  trajectories  by 


x,-tau  I  2  tar.  I  V0  y'>  v')  “  *C'(0 1 

My  [t,y  ,y  )  =  -uTy  (t, y  ,y  )  -  —  exp  — - - - - -J, 


with  the  initial  conditions 


y*a(*o.y'.y')  -v 
i/ta{to,y',y )  =y', 


and  an  energy  by 


Elsa{y',y)  -  y  ([y(a(<,y',y  )]2  +  u2\yta[t,y' ,y  )12)- 
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The  semiclassical  energy  shift  and  uncertainty  for  the  surface  in  the  trajectory  approxima¬ 
tion  are  then 

1  F.ta  „  i  & Eta 

(48) 


Artta  1  d2El\  2  ld2Etg.  .2 
AEs  ~  2~d^~\fc<y  >0+2~d^~^fc<y  >0’ 


and 


(49) 


From  energy  conservation,  the  energy  shift,  A E\*d,  and  energy  uncertainty,  6^E\“a,  for 


the  incidon  are  inferred  to  be 


A Efa  =E‘s  -  A E‘s° 

62E%d 

where  Eg  is  the  initial  mean  energy  of  the  surface.  Using  equations  (42)  and  (43),  the 
momentum  shift  and  uncertainty  for  the  incidon  may  be  written  as 


'ind’ 


(50) 


a  _<a  _  m  a  T7 da  m  r2  nta 

Ap  =  — A Eind  -  -^-6  Eind 


7c 


(51) 


Pfc 
_2 

c2 _<a  ...  m  c2  rpta 
5  P  Jr6  Eind> 

Pfc 

where  p yc  is  the  final  value  of  the  incidon’s  classical  momentum. 

Without  loss  of  generality,  we  may  choose  units  so  that  M  =  w  =  1  and  so  pg/m  =  —1. 
We  may  also  set  Vg  =  1,  since  a  change  in  Vg  is  equivalent  to  a  shift  in  the  coordinate 
x.  Assuming  the  surface  oscillator  is  initially  at  zero  temperature,  the  system  then  has, 
excluding  h,  two  independent  parameters  which  we  take  asf  =  1  j  R  and  rj  =  m/R. 


In  Figures  1-3,  the  exact  momentum  shift  and  that  predicted  by  the  trajectory  ap¬ 
proximation  are  plotted  as  a  function  of  q  for  £  equal  to  1,  2,  and  5,  and  similarly  Figures 
4-6  show  the  momentum  uncertainty  versus  q.  The  ratio  of  the  momentum  uncertainty  in 
trajectory  approximation  to  the  exact  value  is  plotted  in  Figure  7. 


We  see  from  the  figures  that  in  general  the  predictions  of  trajectory  approximation 
deviate  significantly  from  the  exact  result  for  both  the  shift  and  the  uncertainty.  How¬ 
ever,  when  either  r?  or  f  is  small  the  trajectory  approximation  accurately  predicts  the 


momentum  uncertainty  (but  not  the  shift).  In  Appendix  B,  we  prove  that  the  trajectory 
approximation  becomes  exact  fo;  -he  uncertainty  either  if  tj  <  1  or  if  f  <  1  and  that 
t)  <C  1  corresponds  to  the  maximum  deflection  of  the  oscillator  being  small  compared  to 
the  interaction  range,  R,  while  f  <  1  corresponds  to  the  interaction  between  the  incidon 
and  the  oscillator  being  adiabatic.  We  also  note,  without  giving  a  proof,  that  both  the  ex¬ 
act  and  trajectory  approximation  values  for  the  shift  and  uncertainty  have  zeros  whenever 
the  classical  collision  is  perfectly  elastic. 

V.  Conclusions 

Our  main  result  of  equation  (38)  provides  a  simple  method  for  obtaining  the  semiclas- 
sical  limit  of  functions  of  the  incidon ’s  final  momentum.  It  requires  only  the  equilibrium 
correlations  of  the  surface  atoms  and  quantities  derivable  from  classical  trajectories  to  be 
calculated.  We  expect  this  method  to  be  useful  for  physically  interesting  surfaces. 
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equation  (16)  can  be  rewritten 

(d  si  d  ,  d  2. 


/  a  si  a  a  2  a  1  av  u  \  .  . 

\Jt  +  m'a7i  +  +  vS^'ac^5^"’^'0  " 


II(  'C acw)<"']s(fl'‘,c'"')’ 


y/h  a  \ki 


where  the  sum  is  over  all  positive,  odd  G ,  with  G  =  b  where  is  a  nonnegative 

integer;  the  coefficients  J1  are  given  by 


rd-t-iV  ,  a  i  T  OO  r  <i  ,  a  l  -i 

[n  '1 S  ( H  "H*  “ 


where  T  =  YlttJi  with  bu  being  a  nonnegative  integer. 
Using  the  expansion  from  equation  (20) ,  we  obtain 


(f +  mFTi +  iiFri-u‘‘i£l)a^B,‘’Cl‘',)  ~ 1  (ai) 


(I  +  mil  +  ~  “h£-)n(.B„C„ t) 


1  d^ViA^) 


+  2  dA^dA^dA JciBi/Bk dC^B^C^  24dAfldAl/dAjcl 


■  a  A  ^UBv~^9l[BiX'Cli't) 
1  dZyjAj  dZgojB^C^t) 
'  2idAfldAl/dAK]cl  dC^dCudC^ 


The  equation  for  go  is  identical  to  that  for  go,  and  since  g  and  g  have  the  same  initial 


conditions,  we  have  go  =  90-  The  equation  for  g\  is 

/  a  Si  a  a  2  a  \  _  .  _  _  .  c 
lai +  m  dfi  +  <‘*1  ■  1 c“’ 11  - ; 


Subtracting  equation  (A6)  from  (A5),  we  find 
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where  A</i  =  g\  —  g\.  Defining 


Arj(t)  —  J  DBnDCfiAgiiBfijCfi,  t)rj 
Asj(t)  —  J  DBjiDCfx A g\  [Bp ,  C^y  t)sj , 


equation  (A7)  may  be  used  to  show 

A  Sj(t)  =mArj(t) 

assuming  that  g^B^C^)  and  its  derivatives  go  to  zero  as  any  of  its  arguments  goes  to 
infinity.  Since  initially  Arj(£o)  =  Asj(£o)  =  0,  the  unique  solution  to  (A9)  is  A rj(t)  = 
Asj(t)  =  0  which  suffices  to  demonstrate  equation  (26). 

Appendix  B 

Here  we  show  that  the  trajectory  approximation  predicts  the  correct  momentum  un¬ 
certainty  for  either  i|  <  1  or  f  <  1.  We  use  units  with  M  =  u>  =  —p^/m  =  1  and  set 

Vo  =  l. 

Defining  a  set  of  classical  trajectories  by 

Vitr(t,y',y) 

(Bl) 

yir{t,y',y)  =  -  ytr(t,y',y')  - 


with  the  initial  condtions 


xtr(tQ,y',y)  =xQ 

xir{to,y',y)  =  ~  1 
ytr{to,y',y )  -y 
ytr{to,y',y')  =y, 


and  an  energy  by 


Es(t,  y'.y’)  =  5([y'r(e,y’,y')]2  +  [y'r«,y',y  )]2), 


VWIVWlV  v* 


we  have  for  the  exact  semiclassical  energy  uncertainty  of  the  oscillator 

S*ES  =  (^l/c) : 2  <  V2  >0  +  (§7  j/c) 2  <  V1  >0  •  (B4) 

Using  <  y2  >o=<  y2  >0=  V2»  equation  (B4)  may  be  rewritten 

<2  *  - 1{  [(t^)2  -  (f  i/O*]^ +  Kw'"? + (B5, 

From  equation  (49),  an  expression  for  622?£a  can  be  derived  which  is  identical  to  (B5) 
except  that  ytT  is  replaced  by  yia,  while  from  equation  (43)  and  energy  conservation  follow 

«y*  _  £?f£ 

We  first  consider  the  case  with  rj  <  1.  We  may  rewrite  the  second  equation  of  (Bl)  as 
the  integral  equation 

ytr(t )  =  y'  cos(f  -  tQ)  +  y  sin(f  -  *o)  “*  £  /  dt' sin(f  -  t')e^yt  ')~x  ^  )),  {Bl) 

Jto 

which,  using  the  first  equation  of  (Bl),  can  also  be  written 

y*r(t)  =  y'  cos (t  —  to)  +  V  sin(t  —  to)  —  ~  f  dt1  sin(£  —  t')xir{t').  (B8) 

£ Jto 

Taking  y1  and  y1  to  be  infinitesimals,  we  then  have 

ky<r(OI  <  V  f  dt'xtr{t’)  =  r,[xtr{t)  -  xtr{t0)}  <  2v,  { B9 ) 

Jto 

since  |x*r(t)|  <  1.  Therefore,  if  r\  <  1,  then  the  deflection  of  the  oscillator  is  much  less 
than  the  range,  R  = 


If  |?ytr(t)|  <  1,  equations  B(l)  may  be  approximated  by 
yitr(t,y',y') 

y,r{t,y',y’)  =  -  y'r(t, y',y')  - 


(BIO) 


K?' 


Note  that  the  equation  for  xir  is  independent  of  ytr.  Thus 

xtr{t,y',y)  =  xtr{t,  0,0)  =  xd{t),  (Bll) 

and 

S,r(‘,  s' ,  s')  =  -y'r(l,  s',  s')  -  (£12) 

Equation  (B12)  together  with  equation  (45)  imply  that  the  trajectories  use-1  to  calculate 
the  exact  uncertainty  are  the  same  as  those  used  for  the  trajectory  approximation.  Using 
(B7),  we  then  obtain  the  approximate  solution 

ytr{t)  =  yta{t)  =  y'  cos (t  -  to)  +  y'  sin(t  -  to)  -  ?  f  dt'  sin (t  -  t,)e-<**c,(</).  (B13) 

Jto 

Therefore,  from  equation  (B5),  we  have,  for  r\  C  1, 

<2£5  =  =  ||  =  hEtra„„  (B14) 

where  Etrans  is  the  classical  energy  transferred  to  the  oscillator.  Equation  (B14)  together 
with  (B6)  show  that  the  ratio  of  the  trajectory  approximation  momentum  uncertainty 
to  the  exact  momentum  uncertainty  will  be  nearly  unity  if  q  C  1.  It  is  straightforward 
to  verify  that  the  solution  given  in  (B13)  leads  an  energy  shift  of  zero,  and  thus  this 
approximation  tells  nothing  about  the  ratio  of  the  trajectory  approximation  shift  to  the 
exact  shift. 

We  now  treat  the  case  with  f  <.  1.  Defining  \(t)  =  exp[— f**r(t)],  the  second  equation 
of  (Bl)  can  be  rewritten 

ytr(t)  =  -  fA(0e^rtf).  (515) 

The  oscillator’s  motion  will  be  adiabatic  if  A(t)  changes  very  little  over  one  period  of 
oscillation.  Thus  the  adiabatic  condition  is 

j[j||  =  lfi‘r(«)l  «  1-  (£16) 
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Since  |itr(f)|  <  1,  again  considering  y'  and  y'  to  be  infinitesimals,  condition  (B16)  is 
satisfied  if  f  <  1. 

From  the  general  theory  of  adiabatic  motion,14  we  obtain  the  invariant 

HEad,  A)  =  J  [Ead  -  i(a2  +  62)  -  ,  (BIT) 

where 

Ead  =  jlv"'!2  +  |li‘r]2  +  (B18) 

is  the  adiabatic  energy  of  the  oscillator  and  0  is  a  step  function.  In  particular,  I(Eaj,  0)  = 
Ea(}.  Since  both  before  and  after  the  collision  A  is  very  small,  the  invariance  of  I  implies 
that  the  initial  and  final  energies  of  the  oscillator  are  equal. 


In  general,  expanding  Eg  to  second  order  in  y'  and  y'  yields 

Es  =\(yd?  +  |(vc')2 


1 1\(dytr i  ,  ucld2ytr\  i  (dytr i  "i2  i  i>cld2ytr\  L'2 

+5lll7,<v  +  y  ~d^lcl  + t Wlcl)  y  ~w*]cI\y 

l\/dyir.  \2  c,aVr,  (dytr{  \2  ,  .cld2  ytr .  i  .,2 

+  2[(aJ}/'c/)  +y  <?y'2  lc/  +  (  dy'  'c/)  +y  <9y'2  lc/]y 


(519) 


+ 


d2ytr 


\dytr ,  3y*r,  e/  dVr  ,  ,  3y*r,  ,  .e/  cry  , 

.  dy'  ^  ay'  'd  +  y  dy'dy'  'cl  +  dy1  'c/  dy'  'c/  +  y  dy'dy'  'c/ 


y'y 


+  . . . 

Since  in  the  adiabatic  limit,  the  initial  and  final  oscillator  energies  are  equal,  we  have  for 
the  final  energy 


BS  =  iy'2  +  iv'2. 


Comparing  equation  (B19)  with  (B20)  gives 

/V,  \J  /■V.  \2  /V,  'i2.  (3y‘r,  \2 

+  -V57^c)  +(,a7^J  -1 


(£20) 


(521), 


and 


dytr .  ay<r,  ,  3y‘r,  W 


tr 


H/C^rl/c  +  ^rl/c#l/c  =  o. 


<V  ,/c  % 


(£22) 


A  similar  argument  shows  (B21)  and  (B22)  to  be  valid  if  ytr  is  replaced  by  yta.  From 
equation  (B5),  we  thus  obtain  in  the  adiabatic  limit 


62ES  =  62Eisa  =  -(y}c  +  y}c)  =  hEtrans, 


(B  23) 


demonstrating  that  the  trajectory  approximation  gives  the  correct  uncertainty  for  f  <C  1. 
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FIGURE  CAPTIONS 


1)  A p/h  vs.  tj  with  £  =  1  for  the  one-dimensional  example.  The  solid  line  is  the  exact 
value,  and  the  dotted  line  is  the  trajectory  approximation. 

2)  A  p/h  vs.  tj  with  £  =  2.  The  solid  line  is  the  exact  value,  and  the  dotted  line  is  the 
trajectory  approximation. 

3)  A  p/h  vs.  tj  with  £  =  5.  The  solid  line  is  the  exact  value,  and  the  dotted  line  is  the 
trajectory  approximation. 

4)  S2p/h  vs.  rj  with  £  =  1.  The  solid  line  is  the  exact  value,  and  the  dotted  line  is  the 
trajectory  approximation. 

5)  62p/h  vs.  rj  with  £  =  2.  The  solid  line  is  the  exact  value,  and  the  dotted  line  is  the 
trajectory  approximation. 
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6)  S2p/h  vs.  tj  with  £  =  5.  The  solid  line  is  the  exact  value,  and  the  dotted  line  is  the 
trajectory  approximation. 

7)  The  ratio  of  the  momentum  uncertainty  in  the  trajectory  approximation  to  the  exact 
momentum  uncertainty  vs.  rj.  The  solid  line  is  with  £  =  5,  the  dashed  line  is  with  £  =  2, 
and  the  dotted  line  is  with  £  =  1. 
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